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ABSTRACT 


AMRL's  F/  A-18  IFOSTP  fatigue  test  has  introduced  several  new  concepts  in  the 
field  of  aircraft  structural  testing.  One  of  the  innovations  has  been  the 
development  of  the  pneumatic  actuator  or  "rolling  sleeve"  airspring  for  the 
application  of  aircraft  manoeuvre  loads.  These  airsprings  have  been  subjected  to 
several  preliminary  acceptance  test  phases,  one  of  which  simulated  the  effects  of 
high  frequency  buffet  loads.  This  report  details  the  design  and  implementation 
of  the  test  rig  developed  to  conduct  this  type  of  testing. 
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Airsprings 


EXECUTIVE  SUMMARY 


The  RAAF  and  Canadian  Forces  (CF)  have  embarked  on  a  collaborative  full- 
scale  F/A-18  fatigue  test  known  as  the  International  Follow-On  Structural  Test 
Program  (IFOSTP),  with  AMRL  performing  the  testing  of  the  aft  fuselage  and 
empennage,  and  the  Canadians  testing  the  central  fuselage  and  wings. 

The  IFOSTP  project  does  not  use  the  conventional  "whiffle-tree"  method  of 
loading  aircraft  during  full-scale  structural  testing.  Instead,  AMRL  has  designed 
and  developed  another  means  of  applying  aircraft  manoeuvre  loads:-  the 
"rolling  sleeve"  airspring.  These  airsprings  have  a  relatively  low  stiffness  and 
mass  compared  to  a  whiffle-tree  system  and  hence  contribute  less  to  these 
characteristics  of  the  test  article.  The  differential  pressure  in  a  pair  of  rolling 
sleeve  airsprings  distributed  over  the  structure  imparts  the  desired  load 
distribution  into  the  structure.  The  pressurisation  causes  these  inflatable  fabric 
cylinders  to  move  up  and  down  with  a  rolling  action  along  their  support  stands. 

Three  different  phases  of  airspring  testing  have  been  devised:  1)  static  strength 
testing,  2)  durability  testing  under  manoeuvre  loading  and  3)  durability  testing 
under  dynamic  buffet  loading.  This  report  describes  the  development  of  a  test 
rig  for  the  phase  3  testing.  The  test  rig  was  designed  using  a  finite  element 
computer  package  as  well  as  an  analytical  solution  to  confirm  the  results.  The 
working  of  the  analytical  approach  to  this  particular  application  is  one  not 
widely  presented  in  texts  and  so  has  been  described  in  this  report. 
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1.  Introduction 


The  DSTO  Aeronautical  &  Maritime  Research  Laboratory  (AMRL)  has  developed  a 
unique  system  for  the  application  of  manoeuvre  loads  to  aircraft  test  articles.  As  a 
replacement  for  the  more  conventional  hydraulic  jacks  used  in  previous  structural 
tests,  low  stiffness  pneumatic  actuators,  known  as  "rolling  sleeve"  airsprings  (or 
airbags)  have  been  designed  and  developed  to  apply  aircraft  manoeuvre  loads. 
Although  manufacturing  techniques  and  selection  of  suitable  materials  is  an  ongoing 
evolutionary  process,  specific  airspring  units  once  made,  must  satisfy  test 
requirements  to  verify  their  suitability  for  use  in  terms  of  strength,  durability  and 
dynamic  operating  characteristics.  This  airspring  loading  system  will  be  first  trialed 
on  the  empennage  structure  of  the  IFOSTP  full  scale  F/A-18  fatigue  test  article. 
(Reference  1). 

Fatigue  testing  of  the  AMRL  F/A-18  test  article  will  not  only  simulate  manoeuvre 
loading,  but  also  the  high  frequency  dynamic  loading  (or  "buffet")  that  the  aircraft 
experiences  under  certain  flight  conditions.  In  order  to  gauge  how  the  airsprings  will 
react  to  the  differing  load  cases,  independent  tests  have  been  established  to  test 
different  airspring  designs.  Three  tests  have  been  devised;  a  static  strength  test,  a 
manoeuvre  load  simulation  test,  and  a  buffet  simulation  test  which  subjects  the 
airsprings  to  high  frequency  and  small  displacements.  This  report  addresses  the  latter 
requirement. 

Although  AMRL  had  available  a  customised  horizontal  stabilator  rig  that  was 
initially  used  for  developing  the  airspring  loading  system,  for  reasons  of  time,  cost  and 
ease  of  manufacture,  a  simple  cantilevered  beam  structure  was  chosen  for  the 
continuation  of  this  task,  to  be  vibrated  in  its  fundamental  bending  mode  by  an 
electromagnetic  shaker.  It  should  be  noted  that  the  main  aim  in  sizing  the  rig  was  to 
achieve  a  certain  prescribed  tip  deflection  at  the  end  of  the  beam  at  a  prescribed 
frequency.  Another  testing  requirement  of  this  phase  of  test  was  to  determine  the 
damping  characteristics  of  different  kinds  of  airbags,  which  meant  that  the  rig  also  had 
to  cater  for  "free  vibration"  damping  tests. 

This  report  gives  a  brief  background  outlining  the  testing  methodology  developed 
for  AMRL's  F/  A-18  IFOSTP  and  the  need  to  test  the  airbag  system.  It  then  focuses  on 
the  design  of  the  high  frequency,  low  displacement  test  rig,  presenting  the  classical 
approach  for  the  derivation  of  the  equations  of  motion  of  a  cantilever  beam 
undergoing  periodic  motion. 


2.  IFOSTP  Testing  Methodology 


During  F/A-18  flight,  there  are  many  instances  when  significant  empennage 
manoeuvre  loads  occur  simultaneously  with  high  buffet  activity,  causing  structural 
damage  to  the  aircraft's  aft  region  -  see  Reference  2.  In  order  to  accurately  simulate 
this  environment  on  a  test  article,  a  unique  test  system  has  been  developed.  The 
testing  at  AMRL  will  simulate  the  aircraft  manoeuvre  loads  (distributed  aerodynamic 
loads  &  inertia  loads)  by  using  a  new  airbag  (or  airspring)  loading  arrangement 
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instead  of  the  more  conventional  hydraulic  loading  and  contour  board/pad  system 
(collectively  known  as  a  "whiffle-tree"  arrangement).  Because  the  airsprings  have  a 
relatively  low  spring  stiffness  and  mass  compared  to  the  whiffle  tree,  much  less 
stiffness  and  mass  will  be  added  to  the  structure,  minimising  the  effect  on  the  dynamic 
characteristics  of  the  structure,  thus  allowing  the  structure  to  be  excited  near  its 
natural  resonant  frequencies. 

These  airsprings  consist  of  an  opposing  pair  -  one  on  either  side  of  the  structure,  the 
loading  of  the  airspring  pair  being  determined  by  the  differential  pressure  between  the 
airbags.  Air  reservoirs  of  sufficient  volume  are  used  in  conjunction  with  the  airsprings 
to  provide  structural  support  and  to  minimise  the  effective  spring  stiffness  of  the 
system  (stiffness  of  an  airspring  is  inversely  proportional  to  its  volume).  The  dynamic 
buffet  conditions  are  reproduced  using  high  powered,  high  displacement 
electromagnetic  shakers. 

There  are  several  types  of  airsprings  that  will  be  used  for  the  application  of 
manoeuvre  loads.  One  of  these  types  of  actuators  is  the  convoluted  or  "bellows"  type 
commercial  airspring  (see  Figure  1  overleaf)  -  these  are  used  in  the  pneumatic 
suspension  systems  of  trucks  and  for  isolating  machinery  and  equipment  from  other 
structures  in  a  vibrational  environment.  They  are  manufactured  by  Firestone  and 
come  in  various  diameters. 

Development  testing  using  these  commercial  actuators  found  them  inadequate  for 
the  outer  tip  areas  of  the  fins  and  stabilators  as  they  added  too  much  stiffness,  causing 
an  unacceptable  increase  in  modal  frequencies.  The  commercial  airsprings  also  had 
trouble  meeting  the  large  tip  deflections  experienced  during  loading.  This 
necessitated  the  development  of  the  "rolling  sleeve"  airspring/actuator.  The 
development  of  the  rolling  sleeve  airbag  is  detailed  in  Reference  3.  These  low  stiffness 
airsprings  (approx.  1/3  stiffness  of  the  commercial  type)  were  originally  manufactured 
from  a  Kevlar/ MYLAR  sail  cloth  material  and  were  fitted  to  large  volume  reservoirs. 
These  had  very  low  durability  due  to  the  MYLAR  cracking  and  then  abrading  the 
Kevlar. 

Further  development  work  led  to  the  trial  of  different  airbag  materials  in  a  bid  for 
improved  characteristics  such  as  durability,  flexibility,  strength  and  air-tightness.  The 
next  prototype  consisted  of  a  cylinder  of  latex  coated  kevlar  which  was  stitched  and 
sealed.  These  were  followed  by  neoprene  coated  kevlar  that  was  stitched  then  cured 
in  an  autoclave.  Current  airbags  consist  of  a  PVC  coated  polyester  weave  and  have  a 
welded  seam  (high  frequency  heat  bonding  process  that  fuses  the  airbag  coating 
together).  All  bags  are  capped  at  one  end.  Tucks  are  made  around  the  open  end  to 
achieve  a  reduction  in  diameter  to  suit  the  connection  to  the  reservoir  upstand 
support.  The  usual  upstand  support  for  the  airsprings  is  a  length  of  10.79  inch1  (274 
mm)  outer  diameter  steel  tubing  (but  can  vary  depending  on  airbag  diameter). 

Early  experience  from  tests  conducted  at  AMRL  showed  that  failure  of  the  airbags 
tended  to  occur  at  creases  in  the  fabric  caused  by  the  seam  or  tucks  during  the  rolling 
action.  Subsequently,  a  corrugated  plastic  sleeve  was  introduced  to  reduce  the 
creasing  at  the  roll  by  increasing  the  perimeter  of  the  reservoir  upstand  so  that  it  is 
closer  to  that  of  the  airbag.  A  rolling  sleeve  airspring  is  shown  in  Figure  2. 


1 IFOSTP  convention  is  to  use  imperial  units. 
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Figure  1:  Commercial  type  of  airspring  applied  to  horizontal  stabilator. 

The  design  requirements  for  the  airsprings  to  meet  the  IFOSTP  specifications  are: 

Maximum  working  pressure:  40  psi  (275  kPa) 

Peak  System  Pressure  50  psi  (345  kPa) 

Design  maximum  load:  4635  lbf  (20.6  kN) 

Airspring  outer  diameter:  13.5  inches  (343  mm) 

Durability:  120,000  cycles  from  350  lbf  to  1740  lbf  (1 .5  kN  to  7.75  kN) 

2,000,000  cycles  of  1  inch  (25  mm)  stroke  at  1740  lbf  (7.75  kN) 
Stroke:  16  inches  (406  mm) 

Media:  Compressed  air  and  oil  mist. 
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Figure  2:  An  example  of  a  current  rolling  sleeve  airspring. 
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3.  Test  Requirements 


The  airsprings,  when  in  use  during  the  IFOSTP  testing,  will  be  controlled  via  a  control 
law  which  pressurises  the  "loading"  airspring  while  the  opposing  airspring  is  to  be 
kept  at  a  nominal  3  psi  pressure.  During  manoeuvre  loading,  the  deflection  of  the 
flight  surfaces  will  be  in  the  direction  of  loading,  subjecting  the  airsprings  to  a 
combination  of  pressurisation  stressing  and  the  wearing  effects  of  the  rolling  action. 
At  certain  intervals,  while  under  constant  load,  the  structure,  and  hence  airsprings, 
will  be  subjected  to  the  high  frequency  buffet  simulation.  In  these  instances,  "wear”  of 
the  airsprings  will  be  due  to  the  displacement  caused  by  the  rolling  action  only,  and 
not  due  to  major  cyclic  pressure  changes  on  the  body  of  the  airspring. 

Buffet  Loads :  Having  withstood  phase  1  and  phase  2  airspring  development  testing, 
the  airsprings  are  then  subjected  to  the  high  frequency,  small  displacement,  constant 
amplitude  test.  Test  requirements  are  for  airsprings  to  be  inflated  to  apply  a  1740  lbf 
static  load  with  endcap  displacement  of  ±1  inch  for  2,000,000  cycles2- 

To  realise  these  requirements,  the  working  section  of  the  test  rig  was  designed  to 
have  a  natural  frequency  of  approximately  15  Hz  so  as  to  coincide  with  the  15  Hz 
primary  bending  mode  of  the  F/A-18  vertical  tail  -  this  meant  that  minimal  force  input 
would  be  needed  to  drive  the  rig  at  (or  near)  resonance.  The  endcap  displacement  of  ± 
1  inch  was  thought  to  be  sufficient  to  observe  the  effects  of  deflection  on  the  rest  of  the 
airbag.  The  loading  system  used  to  achieve  this  was  an  electromagnetic  shaker  with  a 
1000  lbf  capability. 

Damping:  Another  test  requirement  was  to  measure  the  damping  effects  of  different 
airsprings  to  test  their  suitability  for  use  on  the  IFOSTP  rig.  A  comparative  test  was 
devised  whereby  the  working  section  of  the  test  rig  would  be  made  to  vibrate 
naturally  with  and  without  various  airspring  pairs  attached  -  the  method  of 
determining  the  damping  is  described  in  section  6.2. 

An  airspring  passing  the  above  criteria3  should  be  able  to  withstand  at  least  one  test 
block  of  loading  on  the  IFOSTP  rig. 


2  The  2,000,000  cycles  of  phase  3  testing  was  derived  by  determining  from  F/A-18  operational  data 
(SPEC6G),  the  total  time  spent  when  the  angle  of  attack  (AOA)  >  28°  and  dynamic  pressure  (Q)  >  75  psf, 
and  choosing  rolling  sleeve  airspring  location  47  (stbd)  or  57  (port)  on  the  vertical  tails  as  having  the 
severest  average  displacement  environment  resulting  primarily  from  mode  1  response  to  buffet. 

The  total  time  spent  within  this  AOA-Q  region  was  2828  seconds  per  flight  block.  Since  the  average 
mode  1  frequency  expected  for  the  fatigue  test  is  17  Hz,  and  the  number  of  blocks  for  two  aircraft 
lifetimes  is  40,  the  total  number  of  mode  1  cycles  =  1,923,108.  This  was  rounded  upwards  to  give 
2,000,000  cycles.  From  the  Canadian  Aerospace  Engineering  Test  Establishment  (AETE)  PD88/12  flight 
response  measurements,  the  average  mode  1  response  at  McAIR  measurand  locations  KS16/KT16 
(vertical  tail,  rear  tip,  port  &  stbd  respectively)  for  the  above  AOA-Q  region  was  determined  as  23.76  g 
rms,  and  the  corresponding  response  for  KS01/KT01  (vertical  tail,  forward  tip,  port  &  stbd  respectively) 
was  13.70  g  rms.  Knowing  the  position  of  the  airspring  relative  to  these  transducer  locations,  the  average 
mode  1  response  at  the  airspring  was  determined  to  be  approximately  19  g  rms.  Assuming  this  response 
to  be  a  very  narrow  band  centred  at  17  Hz,  the  expected  peak  displacement  would  be  ±  0.91  inches. 
Hence  an  endcap  displacement  of  ±  1.0  inch  was  chosen  for  this  test. 

3  The  airbags  had  to  survive  at  least  one  full  block  of  loading  (s  300  hrs),  but  preferably  much  longer. 
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Besides  meeting  the  above  criteria,  it  should  be  noted  that  the  construction, 
operation,  performance  and  versatility  of  the  test  rig  were  also  major  points  for 
consideration.  Low  cost  of  new  materials,  availability  of  in-house  equipment  and  ease 
of  manufacture  were  the  driving  factors  behind  the  construction  of  the  rig.  An 
operating  system  that  required  minimal  human  assistance  was  implemented  to  allow 
for  extended  usage  for  operator  free  running  periods  such  as  overnight  running.  The 
test  rig  itself  had  to  perform  predictably  as  well  -  therefore  it  needed  a  solid 
foundation  and  a  working  section  that  was  not  prone  to  erratic  behaviour  or  failure. 


4.  High  Frequency  Test  Rig  Design 


A  cantilever  beam  was  chosen  as  the  design  best  suited  to  fulfil  the  test  requirements. 
However,  as  only  the  fundamental  bending  mode  (in  the  vertical  direction)  was  being 
investigated  in  the  actual  test,  it  was  important  that  the  beam  be  stable  against 
torsional  and  sideways  modes  which  might  occur  at  frequencies  close  to  the  vertical 
mode  (~  15  Hz).  Therefore,  a  rectangular  tube  section  was  used  as  its  characteristics 
best  suited  these  requirements.  A  simple  uniform  cantilever  beam  (no  tip  masses), 
utilising  a  particular  beam  section  that  was  readily  available,  was  chosen  for  the  initial 
design  stages  of  the  rig,  and  its  details  are  shown  in  Figure  3.  The  beam  length  of  70 
inches  is  purely  arbitrary. 


Figure  3  :  Beam  characteristics  used  for  initial  stages  of  design  work. 

Design  of  the  airspring  test  rig  was  performed  using  two  techniques.  One  method 
utilised  the  dynamics  capability  of  the  PAFEC  finite  element  model  (FEM)  package 
and  the  other  method  utilised  an  analytical  solution  approach.  In  order  to  optimise 
the  dimensions  of  the  rig  (to  obtain  the  prescribed  displacements  and  check  stresses), 
sizing  of  the  cantilever  was  performed  using  the  FEM,  as  it  provided  a  quick  and 
straightforward  means  of  experimenting  with  such  parameters  as  beam  geometries, 
material  properties  and  force  inputs. 
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However,  in  the  initial  stages  of  design  work,  an  analytical  solution  for  finding  the 
displacements  of  beams  undergoing  forced  harmonic  motion  was  also  utilised,  so  that 
some  confidence  could  be  gained  from  the  results  given  by  the  FEM.  The  next  two 
sections  present  details  on  both  of  the  methods  used  to  design  the  rig.  A  comparison 
of  the  results  from  the  two  techniques  using  the  beam  described  above  is  presented  in 
section  4.3. 


4.1  PAFEC  FEM  Approach 


The  PAFEC  suite  of  software  was  used  to  solve  for  the  beam  displacements  and 
natural  frequencies.  The  input  file  for  the  FEM  solution  for  an  initial  example  test  run 
is  shown  in  Appendix  1. 

The  geometry  of  the  cantilever  beam  was  defined  along  with  the  type  of  element 
chosen  to  represent  the  structure  in  the  "NODES"  and  "ELEMENTS"  modules 
respectively.  In  this  case  simple  beam  elements  totalling  a  length  of  70  inches  were 
used.  It  should  be  noted  that  the  lengths  and  number  of  these  beam  elements  is 
variable.  If  a  different  beam  section  is  required,  the  dimensions  of  the  beam  elements 
can  be  modified  in  the  "BEAMS"  module  to  represent  different  beam  geometries.  The 
material  properties  of  mild  steel  were  specified  in  the  "MATERIAL"  section.  In  the 
remaining  modules,  the  direction  of  displacement  was  restricted  to  the  vertical  only, 
with  the  built-in  end  of  the  cantilever  being  restrained  in  every  direction.  Also  the 
value  and  location  of  the  applied  force,  the  amount  of  critical  damping  and  output 
information  was  specified  for  various  runs. 

Because  the  time  to  run  a  complete  analysis  on  the  job  was  short  (in  the  order  of  a 
few  minutes),  many  changes  were  made  to  one  parameter  at  a  time,  checking  the 
resulting  displacements,  frequencies  and  stresses. 

For  this  example,  the  force  input  (ie:  position  of  electromagnetic  shaker)  was  set  at  36 
inches  from  the  built-in  end,  corresponding  to  Node  #10  in  the  FE  model.  The 
displacement  output  files  shown  in  Appendix  2  are  for  applied  forces  of  600  lbf  and 
1000  lbf. 


4.2  Analytical  Approach 

This  section  presents  the  analytical  formulation  for  the  vertical  displacements  of  a 
cantilever  beam  undergoing  forced  vertical  vibration,  excited  by  a  concentrated  load 
acting  in  the  vertical  direction  at  an  arbitrary  point  within  the  beam  span.  The 
reference  axes  and  configuration  are  shown  in  Figure  4. 
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The  formulation  of  the  problem  is  presented  below: 

The  differential  equation  (DE)  for  the  flexural  vibration  of  a  beam  has  the  following 
form  (from  Reference  4)  and  is  valid  for  0  <  x  <  1: 


11 

dx2 


3  2y(*,0 

dx2 


+  f(x,t)  =  m(x ) 


d2y(x,t ) 
dt2 


(l) 


with  f(x,t)  being  the  force  per  unit  length.  The  solution  to  the  DE  in  terms  of 
displacement,  y,  as  a  function  of  distance  along  the  beam  (x)  and  time  (t),  is  assumed 
to  be  of  the  following  form: 

y(x,t)  =  Y(x)-F(t)  (2) 

where  Y(x)  represents  the  mode  shape  of  the  beam  and  depends  on  the  spatial 
variable  x  alone,  and  where  F(t)  is  the  generalised  coordinate  indicating  the 
amplitude  of  motion  the  beam  executes  with  time  and  depends  on  t  alone. 

Then  for  the  case  when  /  (x,  t)  =  0  we  have: 


El  Y"" 
m  Y 


(3) 


where  m  is  the  mass  of  the  system  and  go  is  the  frequency  of  F(f),  assuming  that  F(f)  is 
sinusoidal.  Thus,  the  eigenvalue  problem,  associated  with  a  beam  in  bending 
vibration  is  of  the  form  of  the  following  differential  equation: 


dx 2 


EI(x ) 


d2y(x,t) 
dx 2 


=  G)2m(x)Y(x) 


(4) 
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EI(x)  is  referred  to  as  the  flexural  rigidity  and  m(x)  the  mass  per  unit  length  at  any 
point  x  along  the  beam.  In  the  case  of  a  uniform  beam,  El  and  m  are  constant,  and 
equation  (4)  reduces  to: 


d*Y(x) 


rr(jc)  =  o 


where  the  term  k  is  defined  by: 


and  the  constant  a  is  defined  as: 


The  general  expression  for  a  normal  function  can  be  written  as  (from  Reference  5): 

Y;(x)  =  (cos ktx  +  cosh ktx)  +  C2 (cos ktx  -  cosh  ktx) 

+  C3  (sin  ktx  4-  sinh  ktx)  +  C4  (sin  ktx  -  sinh  ktx) 

where  Cv  C2,  C3  and  C4  are  constants  of  integration.  Hence  the  derivatives  of 
I^-(:r)can  be  determined  and  these  are  presented  in  Appendix  3.  The  boundary 
conditions  are: 

(Y)x-_  o  =0 

(n-  =  °  (9) 

cru  =o 

on*  =o 

Using  the  expressions  for  normal  functions  from  Appendix  3  (derived  from  eqn.  (8)) 
and  properties  of  hyperbolic  functions  given  in  Appendix  4,  the  following  frequency 
equation  can  be  obtained  from  (8): 

cos  kl  cosh  kl  =  - 1  (10) 


The  first  six  roots  of  equation  (10)  are  as  follows: 

kj  =  1.875 
k2l  =  4.694 
k3l  =  7.855 
kAl  =  10.996 
k5l  =  14.137 
kj  =  17.279 
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(NOTE:  Should  the  roots  for  higher  modes  be  required,  they  may  be  approximated 
by  the  following  formula  from  the  curve  derived  from  equation  (10)  as  shown  in 
Appendix  5): 


DSTO-TN-OOOl 


k:l  ~  — —  7C  (12) 

2 

Applying  the  boundary  conditions  of  (9),  the  equation  for  the  cantilever  (8)  becomes: 
Yt(x)  =  Bt  [(cos  ktx  -  cosh  kix)  +  ai  (sin  £.x-sinh  kL: c)]  (13) 


where: 


a,  =  — 


cos  &,/+ cosh  kj 
sin  k;l  +  sinh  k:l 


and  constants  Bi  are  selected  to  satisfy  the  conditions  of  normalisation: 


J'  Y*(x)dx  =  1 


ie:  B:  = 


■Jjo  [(cos  k;X  -  cosh  ktx)  +  a,-  (sin  ktx  -  sinh  ktx)f  dx 


It  is  possible  to  obtain  the  expression  for  the  response  y(x,t)  caused  by  the 
concentrated  force  P1  (t),  in  the  form  (See  Reference  5): 

“  Y.Y  rt  P  (t,'\ 

y(x,t)  =  Y-^[  -^sinco ,.(f-r')^'  (17) 

ti  Jo  m 

where  CO,  is  the  i*  natural  frequency,  and  Yn  is  Yt  evaluated  at  the  location  of  Pl . 
Assuming  that  the  vibrations  are  produced  by  a  harmonic  force  P1  (r)  =  P  sin  Qt , 
where  Q  is  the  driving  frequency,  then: 

y(x,  t)  =-£  ]£.=l  ^4[(cos  ktx  -  cosh  ktx)+a  .  (sin  ktx  -  sinh  £,x)] 

x  [(cos  ktxx  -  cosh  ktxx ) + a ;  (sin  fc.jq  -  sinh  kixl  )]Jo  sin  Q.t '  sin  co  t.  (t  -  t')dt' 

~  in  S,=i ^[(cos Kx  ~  cosh ktx)  +  a .  (sin  ktx  -  sinh  /:,.*)] 
x  [(cos  k^  -  cosh  kfa ) + a (sin  ktxx  -  sinh  kix1)^nQ.t-~^n(oit^i 

(18) 

where  the  magnification  factor  (3,  is  given  by: 


l-Q2/©/2 


The  first  part  of  equation  (18)  represents  steady-state  forced  vibrations  of  the  beam, 
whereas  the  second  part  consists  of  transient  free  vibrations.  For  this  particular 
application,  only  the  steady-state  response  will  be  of  practical  importance. 
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Alternatively,  using  ma2  =  El ,  where  a  was  defined  in  (7),  and  the  natural  frequency 
being  given  by  equation  (6),  then: 


(20) 


and  the  response  may  then  also  be  expressed  thus: 

P  °°  B  2l‘i  -| 

y(jt,  t)  =  ~Yj  ^  ^4  [(C0S  kiX  ~  C0Sh  +  a «'  (Sil1  k‘X  ~  Sillh  Mj] 

x  [(cos  ktxx  -  cosh  ktxx  )+af  (sin  kixl  -  sinh  k.xx  )]^sin  fir  -  ^  sin  co  fr  j(3 

(21) 

So,  with  Bif  a i  and  kt  being  the  only  quantities  to  calculate  separately  (the  rest 
being  specific  inputs),  the  displacement  y  at  position  x  and  time  f  can  now  be 
determined.  A  tabulated  form  for  the  solution  of  equation  (21)  is  set  out  in  Appendix 
6. 


In  order  to  calculate  the  displacements  of  the  system,  the  natural  frequencies  also 
had  to  be  determined.  Using  a  variation  of  equation  (20)  above,  the  results  are  listed 
in  Column  8  on  the  table  in  Appendix  6.  The  mass  per  unit  length  m  has  been 
substituted  for  pA/g  and  the  final  result  is  multiplied  by  l/2n  to  give  units  in  hertz. 
The  equation  for  the  first  mode  is  shown  below: 


352  [eJTJ_ 
1  /2  i  pA  '  2k 


(Hz) 


(22) 


It  should  be  noted  that  the  results  from  the  analytical  solution  derived  above  do  not 
account  for  damping  (note  that  the  natural  frequency  of  the  system  is  independent  of 
damping)  -  therefore,  to  keep  things  consistent,  the  FEM  file  (Appendix  1)  also  had 
zero  damping.  The  results  for  two  force  inputs  were  tested  -  600  lbf  and  1000  lbf. 


4.3  Results 


A  comparison  of  tip  deflections  for  two  different  applied  loads  (600  lbf  and  1000  lbf)  at 
several  forcing  frequencies  is  presented  in  Table  1.  The  first  entry  in  each  cell  is  from 
the  FE  solution  and  has  been  taken  directly  from  the  highlighted  sections  of  Appendix 
2.  The  second  entry  (in  parenthesis)  is  from  the  analytical  solution  of  equation  (21) 
and  has  been  calculated  as  per  the  example  in  Appendix  6. 
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Table  1:  Comparison  of  displacements  (in  inches)  from  theFE  and  analytical  solutions. 


£2  =  3  Hz 

£1  =  10  Hz 

£1  =  15  Hz 

£1  =  23  Hz 

£1  =  27  Hz 

£2  =  30  Hz 

£2  =  35  Hz 

P  =  600  lbf 

0.1299 

(0.1299) 

0.1473 

(0.1482) 

0.1800 

(0.1834) 

0.3797 

(0.4171) 

1.299 

(2.2141) 

1.270 

(-0.8532) 

0.2735 

(-0.2404) 

P  =  1000  lbf 

0.2165 

(0.2165) 

0.2455 

(0.2470) 

0.3000 

(0.3056) 

0.6329 

(0.6951) 

2.165 

(3.6899) 

2.117 

(-1.4219) 

0.4559 

(-0.4006) 

One  can  see  that  the  comparison  appears  favourable,  especially  with  values 
extracted  at  frequencies  distant  from  resonance.  For  example,  with  a  shaker  force  of 
1000  lbf  and  a  driving  frequency  of  10  Hz,  the  FE  solution  gives  a  tip  displacement  (at 
node  #19  -  see  schematic  in  Appendix  1)  of  0.2455  inches  while  the  analytical  method 
gives  0.2470  inches.  Also,  some  of  the  analytical  results  are  negative  -  however,  it  is 
the  magnitude  which  is  of  interest. 

Another  point  to  note  is  that  the  analytical  calculation  (eqn.  (22))  produced  a 
fundamental  natural  frequency  of  27.86  Hz  while  the  FEM  solution  calculated  this 
quantity  as  28.51  Hz.  The  FEM  result  of  the  natural  frequency  calculation  is  displayed 
in  the  PAFEC  frequencies  and  displacements  output  file  -  however  this  output  file  is 
very  large  and  has  been  edited  to  only  show  the  displacements  as  presented  in 
Appendix  2.  The  biggest  differences  between  the  two  methods  occurs  near  resonance  - 
however,  regardless  of  which  method  is  used,  the  largest  displacements  are  found  to 
occur  at  this  point  (resonance).  Some  possibilities  that  may  contribute  to  the 
discrepancies  (ie:  differences  in  natural  frequency)  are: 

-  differences  due  to  rounding  errors, 

-  the  accuracy  of  the  FEM  solution  is  affected  directly  by  the  complexity  ofthe 
model.  Errors  are  introduced  depending  on  the  level  of  model  refinement,  as  well 
as  the  choice  of  other  output-analysis  modules. 

With  both  sets  of  data  showing  general  agreement  for  the  simplified  test  case,  this 
indicated  that  the  FEM  results  could  be  used  with  confidence  to  proceed  with  a  more 
detailed  rig  design. 


5.  Test  Rig  Configuration 


The  "final"  version  of  the  cantilever  was  chosen  to  have  a  natural  frequency  of  around 
15  Hz  -  to  coincide  with  the  15  Hz  primary  bending  mode  excitation  of  the  F/A-18 
vertical  tail.  The  provision  for  a  surface  on  which  to  react  the  airspring  pair  was 
needed,  so  it  was  decided  to  make  this  surface  from  half  inch  steel  plate  and  attach  it 
to  the  free  end  of  the  cantilever.  The  dimensions  of  this  plate  were  fixed  -  it  had  a 
certain  thickness  and  it  had  to  be  big  enough  to  support  the  endcaps  of  the  airsprings. 
Consequently,  this  was  a  fixed  parameter  in  the  FEM  (the  plate  was  still  modelled  as  a 
beam  element  but  with  sectional  properties  representative  of  the  plate).  Then  the 
beam  length  was  altered  until  the  target  natural  frequency  of  15  Hz  was  achieved.  The 
same  beam  section  as  described  in  the  example  case  was  used  as  it  provided  good 
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results  and  was  readily  obtainable.  The  position  of  the  electromagnetic  shaker  (ie: 
force  input)  was  also  varied  to  find  an  optimum  location  giving  the  required  tip 
deflections.  This  resulted  in  an  overall  cantilever  length  of  84  inches  -  a  70  inch  beam 
with  a  14  inch  plate  at  the  end.  The  shaker  position  was  at  36  inches  from  the  built-in 
end. 

When  constructing  the  rig,  the  following  method  of  fixation  to  simulate  a 
cantilevered  structure  was  used:  A  2  ton  concrete  block  (48  inches  long,  30  inches 
wide  and  36  inches  high)  was  chosen  as  a  base  on  which  to  attach  the  beam  for  several 
reasons  -  it  provided  a  solid  and  compact  foundation,  as  well  as  being  a  suitable  height 
that  required  little  modification  to  accommodate  the  shaker.  On  top  of  the  concrete 
block  was  bolted  a  one  and  one  quarter  (1.25)  inch  thick  steel  plate.  As  a  result,  when 
it  came  to  ordering  the  beam  from  the  supplier,  an  extra  48  inches  was  added  to  the 
length  of  the  beam,  giving  an  overall  length  of  118  inches.  This  extra  48  inch  section 
was  attached  along  the  whole  length  of  the  upper  surface  of  the  steel  plate.  Pieces  of 
angle-section  were  then  welded  either  side  of  the  beam  and  onto  the  steel  plate.  The 
concrete  block  was  also  anchored  to  the  floor  using  heavy  channel-sections  to  prevent 
motion  induced  during  dynamic  loading. 

The  steel  "pad"  at  the  free  end  of  the  beam  was  actually  20  inches  in  length  -  a  6  inch 
cutout  at  the  end  of  the  beam  was  made  so  that  the  pad  could  be  slotted  in  and  welded 
into  place,  hence  leaving  a  14  inch  overhang. 

A  small  bracket  was  also  welded  to  the  underside  of  the  beam  between  the  built-in 
end  and  shaker.  This  served  as  the  point  from  which  the  beam  could  be  "jacked" 
down  and  released  for  the  free  vibration  damping  tests.  Refer  to  Figure  5. 

Although  not  part  of  the  rig  design,  a  support  frame  to  house  the  airbag  reservoirs, 
mounts,  gauges  and  tubing  was  also  constructed  at  the  free  end  of  the  beam. 


5.1  Test  Rig  Equipment  &  Instrumentation 


The  initial  test  rig  set-up  and  associated  control  devices  consisted  of  the  following 
items: 

TEST  RIG  COMPONENTS 

*  mild  steel  box-beam  section  (6  x  3  x  .236  in,  118  inches  long), 

*  steel  plate  ("end  pad")  -  for  airsprings  to  support  against  (20  x  14  x  0.5  in.), 

*  2  ton  concrete  block  -  as  an  attachment  point  for  cantilever, 

*  steel  base  plate  -  beam  attachment  surface  on  top  of  concrete  block  (48  x  30  x 
1.25  in.), 

*  upright  support  stands  -  for  airspring  and  reservoir  installation, 

*  1  pair  of  rolling  sleeve  airsprings  and  associated  air  reservoirs  &  fixtures, 

*  pressure  transducers  and  gauges  -  to  regulate  and  monitor  pressure  in  airsprings, 

*  hand  operated  hydraulic  jack  and  Canberra  bomb  release  -  to  draw  beam 
downwards  and  release  it  for  free  vibration  damping  tests. 
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Figure  5:  Schematic  of  original  test  rig  set-up. 

LOADING  SYSTEM  &  CONTROLLER 

*  OZY-DYN  signal  generator  -  sets  the  driving  frequency  &  displacement  input 
signal, 

*  LING  1000  lbf  electromagnetic  shaker  &  amplifier  -  the  amplifier  magnifies  the  low 
power  input  from  the  signal  generator  to  an  operational  level  for  the  shaker. 

INSTRUMENTATION 

*  Kistler  accelerometer  -  to  measure  tip  acceleration  with  respect  to  time, 

*  Nicolet  cathode  ray  oscilloscope  (CRO)  -  to  capture  accelerometer  trace, 

*  Wavetek  signal  filter  -  used  to  filter  out  higher  frequency  modes  in  damping  test, 

*  displacement  recorder/plotter  -  provides  a  visual  indication  that  specified 
amplitude  is  being  maintained, 

*  digital  counter  -  to  determine  completion  of  2  million  cycles. 


6.  Testing  Procedure 


6.1  Durability  Testing 

The  pair  of  rolling  sleeve  airsprings  under  test  were  attached  as  shown  in  the  above 
figure  and  pressurised  to  5  psi  before  exciting  the  beam.  For  the  high  frequency 
testing,  the  shaker  "sting"  (the  shaft  transmitting  the  force  output),  was  attached  to  the 
underside  of  the  beam.  The  LING  amplifier  was  then  powered  up  which  in  turn 
started  the  electromagnetic  shaker.  Using  a  frequency  set  control  on  the  signal 
generator  (which  provided  the  inputs  for  the  amplifier),  the  shaker  was  then  set  to 
vibrate  the  beam  near  resonance  (~  15  Hz).  An  accelerometer  positioned  at  the  end  of 
the  beam  (just  in  front  of  the  upper  airspring)  was  used  in  a  feedback  loop  to  prevent 
drifting  of  the  set  frequency.  Another  control  on  the  signal  generator  was  used  to 
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achieve  the  desired  displacements.  Once  the  test  was  running,  this  arrangement 
allowed  for  unmanned  operation,  although  periodic  inspections  were  conducted  to 
ensure  the  equipment  was  operating  normally.  The  rig  showed  no  signs  of  instability. 

The  testing  spectra  consisted  of  cycling  the  airsprings,  with  static  tests  and  load 
versus  pressure  tests  being  performed  at  various  intervals  until  the  2  million  cycles 
were  achieved. 

6.2  Damping  Testing 

When  the  rig  was  used  for  the  damping  tests,  the  shaker  sting  was  disconnected  and 
the  hand  operated  hydraulic  jack  (with  Canberra  bomb  release  attached  to  end)  was 
then  linked  up  to  a  connecting  pin  on  the  underside  of  the  beam.  The  jack  was  then 
used  to  draw  the  beam  downwards  1  inch  before  setting  off  the  bomb  release,  causing 
free  vibration  of  the  beam.  The  displacement  trace  was  captured  on  the  CRO,  using 
the  same  accelerometer  as  mentioned  before.  Measurements  were  taken  of  the  trace  to 
determine  damping  and  frequency.  Damping  was  determined  by  using  the 
displacement  amplitude  ratios  and  chart  shown  in  Figures  6  and  7,  and  as  described  in 
Reference  6.  (For  example,  if  X3  /  xo  =  0.7,  then  read  off  the  n=3  curve  to  obtain  0.019 
or  1.9%  damping). 


Figure  6:  CRO  displacement  trace. 


0.001  0.01  0-! 
FRACTION  OF  CRITICAL  DAMPING 


Figure  7:  Amplitude  ratio  vs.  Fraction  of  critical  damping. 
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The  natural  frequency  of  the  system  was  calculated  by  measuring  the  duration  of 
one  cycle  (in  seconds)  from  the  CRO  display,  then  taking  its  inverse.  (For  example, 
when  the  beam  was  vibrated  without  any  airbags  attached,  there  was  a  time  interval 
of  0.075  seconds  between  xi  and  X2  -  this  gave  a  natural  frequency  of  1/0.075  =  13.33 
Hz.). 


7.  High  Frequency  Test  Results 


The  first  pair  of  airsprings  to  successfully  complete  the  2  million  cycles  were 
constructed  from  a  vinyl  coated  polyester  fabric  (weave  direction  90/0)  manufactured 
by  Polymar  (type  number  8556).  However,  these  particular  airsprings  had  not  been 
subjected  to  phase  2  testing  (they  were  subjected  to  phase  1  testing,  as  are  all 
specimens).  They  were  then  statically  tested  again  to  40  psi  and  showed  no  signs  of 
degradation4.  The  values  of  damping  ranged  from  1.8%  of  critical  at  5  psi  to  2%  of 
critical  at  15  psi. 

The  next  set  of  airsprings  tested  were  made  from  the  same  material  but  were  of  a 
heavier  gauge  (type  number  5407).  These  also  survived  the  test  criteria,  yet  will  not 
be  used  for  IFOSTP  as  they  exhibited  too  much  damping  (from  3.7%  of  critical  at  5  psi 
to  4.05%  of  critical  at  15  psi  -  approximately  twice  the  damping  as  the  #8556  airsprings 
for  this  particular  system). 

As  shown  from  these  tests,  the  airbags  have  performed  well  under  the  phase  3 
conditions.  These  initial  trials  of  the  high  frequency  rig  and  airbags  were  very 
encouraging.  Until  the  full  IFOSTP  test  is  running  and  the  actual  loading  condition  at 
each  airbag  can  be  monitored,  further  refinement  of  the  loading  spectrum  for  the 
durability  tests  is  not  possible. 


8.  Final  Rig  Configuration 


Although  the  design  phase  made  every  effort  to  minimise  the  operating  stress  levels 
throughout  the  beam,  fatigue  cracking  initiated  in  the  weld  attaching  the  beam  to  the 
base  plate  -  consequently,  the  beam  was  replaced  several  times,  with  accompanying 
modifications  to  the  beam-to-base  attachment  being  performed  each  time.  Figure  8 
shows  the  final  rig  configuration,  and  it  should  be  noted  that  the  electromagnetic 
shaker  has  been  replaced  with  an  out-of-balance  rotating  mass  unit,  fixed  to  the  top 
surface  of  the  beam  at  the  free  end.  This  excitation  method  is  driven  by  a  speed 
controlled  DC  motor.  A  review  of  test  methods  led  to  this  change,  which  still 
maintained  unmanned  operation  of  the  rig  (potential  damage  and  costs  to  the 
electromagnetic  shaker  were  high  should  the  beam  fracture  while  vibrating),  and  the 
tip  displacements  were  no  longer  limited  by  the  travel  of  the  shaker  sting.  Also  note 


4  Other  airbags  of  the  same  type  (ie:  #8556)  were  subjected  to  phase  2  tests  using  a  greater  working 
pressure  (2.5  psi  to  35  psi  cyclic  range)  than  that  in  the  original  specification  .  The  durability  of  these 
bags  at  these  pressures  have  so  far  proven  inadequate  for  IFOSTP  requirements.  Further  discussion  on  a 
more  representative  overall  test  spectrum  is  still  ongoing. 
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that  the  method  of  fixing  the  beam  to  the  concrete  block  was  changed  from  a  welded 
attachment  to  a  clamped  arrangement  (steel  channels  over  the  top  of  the  beam  bolted 
to  the  steel  base  plate).  This  serves  to  reduce  stress  concentrations  otherwise  induced 
by  the  welds,  as  well  as  allowing  easier  replacement  of  the  beam. 


9»  Conclusion 


The  Aeronautical  &  Maritime  Research  Laboratory  (AMRL)  has  developed  a  loading 
system  to  apply  manoeuvre  loads  to  the  IFOSTP  test  article  utilising  airspring 
actuators.  As  these  airsprings  are  to  be  subjected  to  the  high  frequency  buffet 
simulation,  a  separate  test  rig  was  commissioned  to  test  individual  airspring  durability 
under  high  frequency  excitation  and  also  to  determine  airspring  damping. 

A  cantilever  beam  structure  was  chosen  as  a  design  for  the  rig,  and  two  methods  of 
sizing  the  rig  were  performed  -  a  finite  element  method  and  an  analytical  approach.  A 
comparison  of  results  from  the  two  methods  proved  favourable,  leading  to  a  final 
design  for  the  test  rig.  The  rig  was  shown  to  serve  its  purpose  well  in  simulating  the 
high  frequency  buffet  conditions  for  testing  of  airspring  actuators. 


Figure  8:  Final  test  rig  configuration. 
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Appendix  1  -  PAFEC  Data  File. 


TITLE  CANTILEVER 
C 

CONTROL 
PHASE=1 
PHASE=  4 
PHASE=6 
PHASE=7 
PHASE=9 
CONTROL . END 
C 

NODES 


NODE . NUM3ER  X 

Y 

1 

0 

0 

2 

4 

0 

3 

8 

0 

4 

12 

0 

5 

16 

0 

6 

20 

0 

7 

24 

0 

8 

28 

0 

9 

32 

0 

10 

36 

0 

11 

40 

0 

12 

44 

0 

13 

43 

0 

14 

52 

0 

15 

56 

0 

16 

60 

0 

17 

64 

0 

18 

68 

0 

19 

70 

0 

C 

ELEMENTS 

NUMBER 

ELEMENT. TYPE 

PROPERTIES 

1 

34000 

1 

2 

34000 

1 

3 

34000 

1 

4 

34000 

1 

5 

34000 

1 

6 

34000 

1 

7 

34000 

1 

8 

34000 

1 

9 

34000 

1 

10 

34000 

1 

11 

34000 

1 

12 

340C0 

1 

13 

34000 

1 

14 

34000 

1 

15 

34000 

1 

16 

34000 

1 

17 

34000 

1 

18 

34000 

1 

C 

BEAMS 

SECTION.  NUMBER  MATERIAL  .NUMBER 

1  12 


TOPOLOGY 

1,2 

2.3 

3.4 

4.5 

5.6 

6.7 

7.8 

3.9 

9.10 
10,11 
11,12 

12.13 

13.14 

14.15 

15.16 

16.17 

17.18 

18.19 


IYY  IZZ 

18.412  6.057 


C 

MATERIAL 

MATERIAL .  NUMBER  E  NU  F.O 

12  29000000  .33  .0007 


C 

RESTRAINTS 

NODE  DIRECTION 


1 

123456 

2 

1345 

3 

1345 

4 

1345 

5 

1345 

6 

1345 

7 

1345 

3 

1345 

9 

1345 

10 

1345 

11 

1345 

12 

1345 

13 

1345 

14 

1345 

15 

1345 

16 

1345 

17 

1345 

18 

1345 

19 

1345 

C 


AREA  ZY  ZZ 

4.025  6.137  4.038 


21 


DSTO-TN-OOOl 


RESPONSE 

TYPE. OF. RESPONSE  OUTPUT. TYPE 

1  1 
C 

DAMPING 

MODE=0 

FREQUENCY  DAMPING . RATIO 

0  0.0 

35  0.0 

C 

FREQUENCIES . FOR . ANALYSIS 
TYPE  START  FINISH  STEP 

1  0  35  1 

C 

SINE. LOADING 

NODE. NUMBER  DIRECTION . OF. LOAD  TABLE. NUMBER 

10  2  1 

C 

TABLE. OF. APPLIED. FORCES 
TABLE  BASIS. VALUE  VALUE. LIST 

1  0  500  0 

1  35  600  0 

C 

SINUSOIDAL .OUTPUT 
NODE  DIRECTION 

10  2 

19  2 

C 

MASTERS 

NODE. NUMBER  DIRECTION 

10  2 

19  2 

C 

MODES .AND . FREQUENCIES 
AUTOMATIC . MASTERS 
SO 
C 

END. OF. DATA 


SCHEMATIC  OF  CANTILEVER  BEAM  MODEL. 
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Appendix  2. 

SECTION  OF PAFEC  DISPLACEMENT  OUTPUT  FILE  (P  =  600  lb/). 


. 

NODE  10.  DIRECT  Z ON  2 

NCOS  19. 

DIRECTION  2 

LOAD . • 
CASE  * 

FREQUENCY 

DISPLACEMENT  * 

PHAS3. ANGLE 

DISPLACEMENT 

PHASE. ANGLE 

1  * 

.000 

.S312S-01  * 

360.0 

.1284 

360.0 

1  * 

1.000 

. 5313E-01  • 

360.0 

.1235 

360.0 

1  * 

2.000 

.5336E-01  * 

360.0 

.1290 

350.0 

1  • 

3 .000 

. 53  56E-01  * 

360.0 

ri2«"i 

360.0 

1  ' 

4 . 000 

- 

3  oQ . 0 

.  i  3 11 

3  60.0 

1  • 

5.000 

.S465E-01  * 

360.0 

.1327 

360.0 

X  * 

6.000 

. 5535E-01  * 

3  60.0 

.1346 

360.0 

1  * 

7.000 

.56223-01  * 

360.0 

.1370 

360.0 

1  * 

3.000 

.57243-01  ' 

360.0 

.1399 

360.0 

1  * 

9.000 

.58463-01  • 

360.0 

.14  3  3 

360.0 

1  • 

10.000 

.59393-01  • 

360.0 

_ Eim.- 

3  60.0 

1  * 

11.000 

"I" 

.0-533-01  * 

360  . U 

.1320 

360 .0 

1  - 

12.000 

.  63  50E-01  " 

360.0 

.1574 

360.0 

1  • 

13.000 

.  6577E-0 1  • 

360.0 

.1638 

360.0 

1  * 

14.000 

.  5343  E-0 1  • 

3  60.0 

.1712  , 

360.0 

i  * 

15.000 

. 7157E-0 1  • 

350.0 

1  .1800  1 

360.0 

1  * 

16.000 

— .25273-ar*- 

JaO 

.1904 

3oO 

1  * 

17.000 

.79703-01  * 

360.0 

.2028 

360.0 

1  ' 

13.000 

.35053-01  * 

3  60.0 

.2178 

360.0 

1  * 

19.000 

.91503-01  * 

350.0 

.2361 

360.0 

1  ’ 

20.000 

.99773-01  * 

3  60.0 

.2591 

360.0 

1  * 

21.000 

.1102 

3  50.0 

.2334 

360.0 

1  • 

22.000 

.1240 

360.0 

.3269 

360.0 

i  * 

23.000 

.1423 

360.0 

1 

360.0 

1  * 

24.  OOF 

.1701  ~ 

jaU  .  0 

.  4do4 

j  60 . 0 

1  * 

25.000 

.2131 

360.0 

.5770 

360.0 

1  * 

26.000 

.2905 

3  60.0 

r-I242_, 

360.0 

1  * 

27.000 

.4703 

360.0 

11.299  | 

360.0 

1  * 

28.000 

T7T53 - - 

IdFS 

3.767 

3  60.0 

1  * 

29.000 

1.397 

130.0 

3 . 942 

180.0 

1  * 

30.000 

.  4450 

130.0 

130.0 

1  * 

31.000 

~  .2593  ' 

130. 0 

.  *  4  9^3 

UdVO 

1  * 

32.000 

.1303 

180.0 

.5276 

180.0 

1  * 

33.000 

.1366 

130.0 

.4049 

180.0 

1  * 

34.000 

.1089 

130.0 

,  0212-, 

180.0 

* 

i  * 

35.000 

.3978E-01  * 

130.0 

1.2735  1 

_ L 

180.0 

SECTION  OFPAFEC  DISPLACEMENT  OUTPUT  FILE  (P  =*  1000  lb/). 


. 

. 

NODE  10.  DIRECTION  2 

NODE  19.  DIRECTION  2 

LOAD.* 
CASS  * 

FREQUENCY  *' 

DISPLACEMENT  * 

PHASE.  ANGLE  * 

DISPLACEMENT 

PHASE  .ANGLE 

1  * 

1  * 

1  * 

1  * 

.000  - 
1.000  * 
2.000  * 
3.000  * 

. 8854E-01  * 

.  38 54 E-0 1  • 

. 3394E-01  * 
.3944S-01  " 

360.0 

360.0 

3  50.0 

360.0 

.2140 

.2142 

.2151 

mm 

360.0 

360.0 

360.0 

360.0 

1  * 

1  * 

1  * 

1  * 

1  * 

1  * 

1  * 

4.000  * 
5.000  * 
6.000  * 
7.000  * 
8.000  * 
9.000  * 
10.000  * 

.90x33-01  " 

. 9109E-Q 1  * 
.92263-01  * 
.93693-01  * 
.95413-01  * 
.9743E-01  * 
.99813-01  * 

360.0 

360.0 

360.0 

360.0 

360.0 

3  60.0 

360.0 

.  2135 
.2211 
.2244 
.2284 
.2332 
.2338 

1  .2455  | 

360 .0 

360.0 

360.0 

360.0 

360.0 

360.0 

360.0 

1  * 

1  * 

1  * 

1  * 

n.ooo-* 
12.000  * 
13.000  * 
14.000  ♦ 
15.000  * 

.102  = 

.1053' 

.1095 
.1141 
.119  3 

J  oO  .  0 

360.0 

360.0 

360.0 

3  60.0 

,233j 
.2624 
.2730 
.2354 
( .3000  1 

360 .0 

360.0 

3  60.0 

3  60.0 
-  360.0 

1  * 
1  * 

1  * 

1  * 
1  * 

1  * 
1  * 
1  • 

16.000  * 
17.000  * 
13.000  * 
19.000  • 
20.000  * 
21.000  * 
22.000  * 
23.000  * 

.1265 

.1328 

.1417 

.1527 

.1563 

.1337 

.2065 

.2330 

UUU  txiu  uwu 

a>  Oi  o»  o»  o» 
oooooooc 

oooooooc. 

■  ii/3 
.3380 
.  3630 
.3936 
.4313 
.4806 
.  5443 

1 . 6329  | 

uuuuuuwu 
Ok  a»  oi  cr»  C7i  cn  o 
OOOOOOOC 

oooooooc 

1  * 
1  * 
1  * 
1  . 

24.000  • 
25.000  * 
25.000  * 
27.000  * 

.  3552 
.4341 
.7339 

3  60.0 

3  60.0 

3  60.0 

.  /  oOo 
.9617 

,1-124., 
i  2 . 165  1 

*  JOU . u 

*  360.0 

*  350.0 

*  360.0 

1  * 
1  * 

1  * 

28.000  * 
29.000  • 
30.000  * 

l.iSO 

2.323 

.7417 

~  360.0 

180.0 

130.0 

6.279 

6^70 
12.117  1 

*  3  oO  .  0 

*  130.0 

*  130.0 

1  * 
1  * 
1  * 
1  * 
i  . 

31.1300  - 
32.000  * 
33.000  * 
34.000  * 
35.000  * 

.  4322 
.3005 

.2277  * 

1 1313 
.1495 

- TJO - 

130.0 

130.0 

130.0 

130.0 

TT3T5 
.3793 
.6749 
«  4  c  1 

1.4559  1 

“  130.0 

*  180.0 

*  130.0 

180.0 

*  130.0 
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Appendix  3. 

Yt  =  C[(cos^x  +  cosh ktx)  +  C2{cosk,x  -  cosh ktx) 

+  C3 (sin  ktx  +  sinh  ktx)  +  C4 (sin  ktx  -  sinh  k,x) 

Y/  =  Cykt{- sin  k.x  +  sinh  k.x)  +  C2k.(- sin  k}x  -  sinh  k.x) 

+  C3^.(cosA:(x  +  cosh&;x)  +  CJct{cosktx  -  cosh ktx) 

Y"  =  Clkf(-coskix  +  coskkix)+C2kf(-coskix-cosh.kix) 
+  C,kf  (-  sin  k.x  +  sinh  ktx)  +  C4kf  (-  sin  k,x  -  sinh  k,x) 

Y.”  -  Cxk]  (sinAix  +  sinh  k.x)  +  C2k]  ( sin  k.x  -  sinh  ktx) 

+  CJc]  (-  cosksx  +  cosh k{x)  +  C4k]  (-  cosk,x  -  cosh ktx) 

Appendix  4. 

(coshx)^  =  [k-kf)  =  1 

j-0 

(smh*)„  =  (^)  =0 

x=0 

coshx2  -sinhx2  =  1 

Appendix  5. 


cosp£-coshf3Z,  =  -l 

cosp  L  =  ^- - =  (f)2 

coshpZ, 

The  above  equation  is  satisfied  atnjq  =  <j>2 


ANAL  YTICAL  SOLUTION  OF  EQUA  TION  (21) 

29  x  106  psi  1  =  6.057  in4  Beam  length, x  -  l  =  70  in  Shaker  @  x,  =  36  in  A  =  4. 025  in2  p  =  0. 283  lb  /  in1  g  =  3 86  in  /  s 
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Appendix  6. 
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Tlie  sum  oftlie  products  (122100. 5193)  is  then  multiplied  by  P/EI  (3.416  x  IOE-6)  to  give  the  displacement  which  is  0.4171  inches. 
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